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1 A triple point 



Through the last fifty years, the possible existence and implications of primordial mag- 
netic fields became a very useful cross-disciplinary area at the interface of cosmology, 
astrophysics and high-energy physics|I|, |^, |^, ^, ^, ^, 0. From astrophysical observations, 
we do know that planets, stars, the interstellar medium and the intergalactic medium are 
all magnetized. The magnetic fields in these environments have values ranging from the 
/X G of the intra-cluster medium, to the G (in the case of the earth) up to presumably 
10^^ G, the typical magnetic fields of neutron stars. In principle, observations of magnetic 
fields in galaxies (and clusters of galaxies) could discriminate between a direct primordial 
origin of the large-scale fields and a primordial origin mediated by a dynamo amplification. 
None of the two options are, at the moment, supported by clear observational evidence. 
Furthermore, in both approaches there are various theoretical assumptions which have 
been (but still need to be) carefully scrutinized. 

In high-energy physics, the possible existence of intergalactic magnetic fields is one 
of the crucial unknowns in the analysis of ultra-high energy cosmic rays above the GZK 
cut-off. The interplay of high-energy physics and astrophysics is indeed present since the 
origin of this subject. In 1949 the scientific argument between Fermi ^j, on one side, and 
Alfven |T0|, Richtmyer and Teller []1T|, on the other, concerned exactly the possible 



existence of galactic magnetic fields. Fermi was convinced that high-energy cosmic rays 
are in equilibrium with the whole galaxy while Alfven was supporting the idea that high 
energy cosmic rays are in equilibrium with stars. In order to make his argument consistent, 
Fermi postulated (rather than demonstrated) the existence of a /i G galactic magnetic 
field. Fermi thought that the origin of this field was primordial. 

In cosmology the possible existence of magnetic fields prior to decoupling can influence 
virtually all the moments in the thermodynamical history of the Universe. Big-bang 
nucleosynthesis (BBN), electroweak phase transition (EWPT), decoupling time are all 
influenced by the existence of magnetic fields at the corresponding epochs. If magnetic 
fields were originated in the past history of the Universe, their birth should be related, in 
some way, to the interplay of gravitational and gauge interactions. Superstring theories 
and higher-dimensional theories, formulated through the past thirty years, pretend to give 
us some hints on the possible form of such an interplay and their implications may be 
useful to consider. 

The present paper is organized as follows. In Section II the basic ideas on large scale 
magnetic field structure and observations will be briefly outlined. Section III collects some 
considerations on the evolution of magnetic fields. In Section IV the problem of the origin 
will be illustrated with particular attention to models where there is an effective evolution 
of the gauge coupling. Section V deals with the possible implications of hypermagnetic 
fields for the EW physics and for the generation of the BAU. In Section VI it will be shown 
that if hypermagnetic fields are present at the EW epoch, matter-antimatter fluctuations 
are likely to be produced at BBN. In Section VII the implications of hypermagnetic fields 
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for the GW backgrounds at the LISA and VIRGO/LIGO frequencies will be discussed. In 
Section VIII some speculations on the possible Faraday rotation of the CMB polarization 
will be presented. Section IX contains some concluding remarks. 



2 A primer on magnetic fields observations 

There already excellent reviews on the measurements of large scale magnetic fields in 
diffuse astrophysical plasmas 0, ||, |, ^. Here we want just to give an elementary 
primer on the main ingredients of the detection strategies and focus the attention on 
recent observations. 

2.1 Local and global observables 

In order to measure large scale magnetic fields, one of the first effects coming to mind, 
is the Zeeman splitting. The energy levels of an hydrogen atom in the background of a 
magnetic field are not degenerate. The presence of a magnetic field produces a well known 
splitting of the spectral lines: 



From the estimate of the splitting, the magnetic field intensity can be deduced. Now, the 
most common element in the interstellar medium is neutral hydrogen, emitting the famous 
21-cm line (corresponding to a frequency of 1420 MHz). Suppose that a magnetic field of 
fi G strength is present in the interstellar medium. From Eq. the induced splitting, 

Auz ~ 3Hz, can be estimated. Hence, Zeeman splitting of the 21-cm line generates 
two oppositely circular polarized spectral lines whose apparent splitting is however sub- 
leading if compared to the Doppler broadening. In fact, the atoms and molecules in the 
interstellar medium suffer thermal motion and the amount of induced Doppler broadening 
is roughly given by 



where Vth is the thermal velocity oc yjT/m where m is the mass of the atom or molecule. 
The amount of Doppler broadening is Az/Dop ~ 30 kHz which is much larger than the 
Zeeman splitting we ought to detect. Zeeman splitting of molecules and recombination 
lines should however be detectable if the magnetic field strength gets larger with the 
density. Indeed in the interstellar medium there are molecules with an unpaired electron 
spin. In these cases a Zeeman splitting comparable with the one of the 21-cm line can be 
foreseen. These molecules include OH, CN, CH and some other. In the past, OH clouds 
were used in order to estimate the magnetic field (see and references therein). The 
possible caveat with this type of estimates is that the measurements can only be very 
local. The above mentioned molecules are much less common than neutral hydrogen and 
are localized in specific regions of the interstellar medium. 




(2.2) 
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The first experimental evidence of the existence of large scale magnetic fields in galaxies 
came from the synchrotron emission. The emissivity formula for the synchrotron depends 
upon B± and upon the relativistic electron density. The synchrotron has an intrinsic 
polarization which can give the orientation of the magnetic field, but not the specific 
sign of the orientation vector. The relativistic electron density is sometimes estimated 
using equipartition, i.e. the idea that magnetic and kinetic energy densities may be, 
after all, comparable. Equipartition is not an experimental evidence, it is a working 
hypothesis which may or may not be realized in the system under observation. For 
instance equipartition probably holds for the Milky Way but it does not seem to be valid 
in the Magellanic Clouds [l^. The average equipartition field strengths in galaxies ranges 



from the 4/i G of M33 up to the 19/i G of NGC2276 [|14 . 

In order to infer the magnitude of the magnetic field strength Faraday effect has been 
widely used. When a polarized radio wave passes through a region of space containing a 
plasma with a magnetic field the polarization plane of the wave gets rotated by an amount 
which is directly proportional to the square of the plasma frequency (and hence to the 
electron density) and to the Larmor frequency (and hence to the magnetic field intensity). 
Calling (j) the shift in the polarization plane of the wave, a linear regression, connecting 
the shift in the polarization plane and the square of the wavelength of observation, can 
be obtained: 

= 00 + RMAl (2.3) 

By measuring this relation for two (or more) separate (but close) wavelengths, the angular 
coefficient of the regression can be obtained and it turns out to be 

in units of rad/m^ when all the quantities of the integrand are measured in the above 
units. As reminded, RM should be performed at sufficiently close wave-lengths. Typically 
the angles should be determined with an accuracy greater than 50 ~ ivr. Otherwise 
ambiguities may arise in the determination of the angular coefficient appearing in the 
linear regression of Eq. ( |2.3| ) [|T], Q . 

It should be appreciated that the RM contains not only the magnetic field (which 
should be observationally estimated), but also the column density of electrons. From the 
radio-astronomical observations, different techniques can be used in order to determine the 
column density of electrons. One possibility is to notice that in the observed Universe there 
are pulsars. Pulsars are astrophysical objects emitting regular pulses of electromagnetic 
radiation with periods ranging from few milliseconds to few seconds. By comparing the 
arrival times of different radio pulses at different radio wave-lengths, it is found that 
signals are slightly delayed as they pass through the interstellar medium exactly because 
electromagnetic waves travel faster in the vacuum than in an ionized medium. Hence, 
from pulsars the column density of electrons can be obtained in the form of the dispersion 
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measure, i.e. DM oc / UeM. Dividing the RM by DM, an estimate of the magnetic field 
can be obtained. 

Aheady this simple-minded account of the main experimental techniques used for 
the detection of large scales magnetic fields shows that there may be problems in the 
determination of magnetic fields right outside galaxies. There magnetic fields are assumed 
to be often of nG strength. However, due to the lack of sources for the determination of 
the column density of electrons, it is hard to turn the assumption into an experimental 
evidence. 

Finally, an interesting source of observational informations on galactic magnetic field 
structure is the polarization of synchrotron emission. Magnetic fields in external galaxies 
have an uniform and a random component. Synchrotron polarization at high radio fre- 
quencies (where Faraday rotation is small) can be used in order to estimate the relative 
weight of the mean and random parts of a given field since the polarization essentially 
depends upon the ratio between the uniform and the total (i.e. random plus uniform) 
magnetic field 0. 

Since various theoretical speculations suggest that also clusters are magnetized, it 
would be interesting to know if regular Abell clusters posess large scale magnetic fields. 
Different results in this direction have been reported [^, |16|, Some studies during 



the past decade IT^, |l6l dealt mainly with the case of a single cluster (more specifically 
the Coma cluster). The idea was to target (with Faraday rotation measurements) radio 
sources inside the cluster. However, it was also soon realized that the study of many radio 
sources inside different clusters may lead to experimental problems due to the sensitivity 
limitations of radio-astronomical facilities. The strategy is currently to study a sample of 
clusters each with one or two bright radio-sources inside. 

In the past it was shown that regular clusters have cores with a detectable component of 



RM |T^, |TB| . Recent results suggest that jj, Gauss magnetic fields are indeed detected inside 



regular clusters |T^. Inside the cluster means in the intra-cluster medium. Therefore, 



these magnetic fields cannot be associated with individual galaxies. 

Regular Abell clusters with strong x-ray emission were studied using a twofold tech- 



nique |]19|, From the ROSAT Q full sky survey, the electron density has been de- 
termined. Faraday RM (for the same set of 16 Abell clusters) has been estimated through 
observations at the VLA The amusing result (confirming previous claims based only on 



one cluster [^, |T6[) is that x-ray bright Abell clusters possess a magnetic field of ^ Gauss 
strength. The clusters have been selected in order to show similar morphological features. 
All the 16 clusters monitored with this technique are at low red-shift {z < 0.1) and at 
high galactic latitude (|6| > 20°). 

These recent developments are rather promising and establish a clear connection be- 

*The ROetgen SATellite was flying from June 1991 to February 1999. ROSAT provided a map of the 
x-ray sky in the range 0.1-2.5 keV. For the ROSAT catalog of X-ray bright Abell clusters see [ p2| . 

^The Very Large Array telescope is a radio-astronomical facility consisting of 27 parabolic antennas 
spread 20km^ in the New Mexico desert. 
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tween radio-astronomical techniques and the improvements in the knowledge of x-ray 
sky. There are various satellite missions mapping the x-ray sky at low energies (ASCA, 
CHANDRA, NEWTON []) . There is the hope that a more precise knowledge of the surface 
brightness of regular clusters will help in the experimental determination of large scale 
magnetic fields between galaxies. 

It is interesting to notice that intra-cluster magnetic fields of G strength can in- 
duce Faraday rotation on CMB polarization. By combining informations from Sunyaev- 
Zeldovich effect and X-ray emission from the same clusters, it has been recently suggested 



that a richer information concerning electron column density can be obtained ||23|. In 
Fig. m the results reported in |jl9| are summarized. In Fig. |l] the RM of the sample of 
x-ray bright Abell clusters is reported after the subtraction of the RM of the galaxy. At 
high galactic latitude (where all the observed clusters are) the galactic contribution is 
rather small and of the order of 9.5rad/m^. In Fig. the open points represent sources 
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Figure 1: From Ref. the RM deduced from a sample of 16 X-ray bright Abell clusters 
is reported as a function of the source impact parameter. 

viewed through the thermal cluster gas, whereas the full points represent control sources 
at impact parameters larger than the cluster gas. The excess in RM attributed to clusters 
sources is clearly visible. 

Using the described techniques large scale magnetic fields can be observed and studied 
in external galaxies, in clusters and also in our own galaxy. While the study of external 
galaxies and clusters may provide a global picture of magnetic fields, the galactic obser- 
vations performed within the Milky Way are more sensitive to the local spatial variations 
of the magnetic field. For this reasons local and global observations are complementary. 



*ASCA is operating between 0.4 AND 10 keV and it is flying since February 1993. CHANDRA (NASA 
mission) and NEWTON (ESA mission) have an energy range comparable with the one of ASCA and were 
launched, almost simultaneously, in 1999. 
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The flipped side of tlie coin, as we will see in the second part of the present Section, is 
that the global structure of the magnetic field of our galaxy is not known directly and to 
high precision but it is deduced from (or corroborated by) the global knowledge of other 
spiral galaxies. 



2.2 Geometrical models 



Since the early seventies Q the magnetic field of the Milky way was shown to be 
parallel to the galactic plane. RM derived from pulsars allow consistent determinations 
of the magnetic field direction and intensity p8| , |29[] . In the Milky Way, the uniform 
component of the magnetic field is thought to lie in the plane of the galactic disk and 
it is thought to be directed approximately along the spiral arm. There is, though, a 
slight difference between the northern and southern hemisphere. While in the southern 
hemisphere the magnetic field is roughly 2/iG, the magnetic field in the northern emisphere 
is three times smaller than the one of the southern hemisphere. Differently from other 
spirals, the Milky Way has also a large radio halo. Finally RM data seem to suggest 
that the magnetic field flips its direction from one spiral arm to the other. As far as the 
stochastic component of the galactic magnetic field is concerned, the situation seems to 



be, according to the reported results, still unclear p5|, |2^, It is, at present, unclear 
if the stochastic component of the galactic magnetic field is much smaller than (or of 
the same order of) the related homogeneous part. The geometries of large scale magnetic 
fields in spiral galaxies can be divided into two classes. We can have axysymmetric spirals 
(ASS) and bisymmetric spirals (BSS). This classification scheme refers, respectively, to 
even and odd parity with respect to rotation by an angle n around the galactic center. 
Speculations concerning the origin of galactic fields prefer to associate a primordial field 
with a BSS configuration while the ASS would be more associated with a field produced 
through a strong dynamo activity. 

In the case of the Milky way, as we saw, the magnetic field flips its direction from 
one spiral arm to the other and then, as pointed out by Sofue and Fujimoto (SF) the 
galactic magnetic field should probably be associated with a BSS model. In the SF model 
the radial and azimuthal components of the magnetic field in a bisymmetric logarithmic 
spiral configuration is given through 

Br = f(r) cos ( 6 — f3ln — ) sinp 
V roJ 

Be = fir)cos(^e- pln^^cosp (2.5) 

where tq ~ 10.5 kpc is the galactocentric distance of the maximum of the field in our 
spiral arm, j3 = 1/ tanp and p = 10° is the pitch angle of the spiral. The smooth profile 
/(r) can be chosen in different ways. A motivated choice is [^, ^ 



/(r) = 3^tanr ( — )/iG (2.6) 
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where ri = 8.5 kpc is the distance of the Sun to the galactic center and r2 = 2 kpc. 
The original model of SF does not have dependence in the z direction, however, the z 
dependence can be included and also more complicated models can be built . Typically, 
along the z axis, magnetic fields are exponentially suppressed as exp \—zIzq\ with zq ^ 4 
kpc. The structure of magnetic fields can be relevant when investigating the propagation 
of high-energy protons Bl, as noticed already long ago [p3] (see also pHl E^, RO]). 



3 Magnetic field evolution (s) 



The galaxy is a gravitationally bound system formed by fluid of charged particles which 
is globally neutral for scales larger than the Debye sphere. In the interstellar medium, 
where the electron density is approximately 3 x 10~^ cm~^, the Debye sphere has a radius 
of roughly 10 m. Moreover, the galaxy is rotating with a typical rotation period of 
3 X 10*^ yrs. Two complementary descriptions of the plasma can then be adopted. The 
first possibility is to study full kinetic system (the Vlasov-Landau equations [^, p9|). 
The second (complementary) description relies on the magnetohydrodynamical (MHD) 
treatment. Among the discussions of MHD effects, dynamo theory is particularly relevant. 
A full account of the various aspects and debates concerning the dynamo theory can be 
found in excellent textbooks Q and reviews After some elements of the Vlasov- 



Landau description, the the basic ideas concerning dynamo theory will be introduced. 



3.1 Elements of a kinetic discussion 



Already in flat space ||4^, and, a fortiori, in curved space |^T|, the kinetic approach 



is important once we deal with electric fields dissipation, charge and current density 
fluctuations and, in more general terms, with all the high frequency and small length 
scale phenomena in the plasma [^, ^ . The few elements of kinetic description will be 
given directly in curved spaces since they may be relevant for some applications [Q. 

Consider a conformally fiat Friedmann-Robertson- Walker (FRW) metric written using 
the conformal time coordinate 



ds'^ = a^{ri)[dri^ - dx% (3.1) 

Furthermore, consider an equilibrium homogeneous and isotropic conducting plasma, 
characterized by a distribution function fo{p) common for both positively and negatively 
charged ultrarelativistic particles (for example, electrons and positrons) . Suppose now 
that this plasma is slightly perturbed, so that the distribution functions are 

/+(f,p,r/) = /o(p) + (5/+(f,p,r/), /_(f,p,r/) = /o(p) + (5/„(f,p,r/), (3.2) 

where + refers to positrons and — to electrons, and p is the conformal momentum. The 
Vlasov equation defining the curved-space evolution of the perturbed distributions can be 
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written as 



"441 



df+ df+ , 5n 



drj 
drj 



+ v 



dx 
df- 

dx 



e{E + vx B) 



df- 



\ drj / coll 

iir) ■ 

V drj / coll 



(3.3) 
(3.4) 



where the two terms appearing at the right hand side of each equation are the colhsion 
terms. The electric and magnetic fields are rescaled by the second power of the scale 
factor. This system of equation represents the curved space extension of the Vlasov- 
Landau approach to plasma fiuctuations All particle number densities here are 

related to the comoving volume. By subtracting Eqs. ( p.3| ) and ( |3.4| ) we obtain the 
equations relating the fiuctuations of the distributions functions of the charged particles 
present in the plasma to the induced gauge field fluctuations: 

d d dfo 

Q:^f{x,P,t) + v- —f{x,p,t) + 2eE ■ = -u{p)f, 



V ■E = 

V X E + B' 

V - 5 = 0, 

V X B-E' 



£pf{x,p,rj), 



0, 



d^pvf{x,p,rj), 



(3.5) 



where / (x, p, fj) 

El- 



6f+{x,p,rj) — 6f-{x,p,rj) and z/(p) is a typical frequency of collisions 



Now, if Sf±{x,p,rj) 7^ at the beginning of the radiation dominated epoch rjo and 
E{x, rjo) ~ B{x, Tjo) = initially, the magnetic field at later times can be found from Eqs. 

Various useful generalizations of the Vlasov-Landau system to curved spaces 



is given in 45, 46, 47 . 



3.2 Effective (MHD) description and dynamo instability 

For scales sufficiently large compared with the Debye sphere and for frequencies sufficiently 
small compared with the plasma frequency, the spectrum of plasma excitations obtained 
from the kinetic theory matches the spectrum obtained from an (effective) MHD descrip- 
tion Furthermore, since the galaxy is rotating and since the conditions of validity of 



the MHD approximation are met, it is possible to use the so-called dynamo instability in 
order amplify a small magnetic inhomogeneity up to the observed value. This is, at least, 
the hope . 

The pioneering attempts towards a MHD description of interstellar plasma, go back 
to the works of Alfven [0] and of Fermi and Chandrasekar [T2\. Since the work of Parker 



48| , on the so-called a — Q theory, the dynamo effect has been used in order to explain 
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(or, to ease) the problem of the origin of galactic magnetic fields. The standard dynamo 
theory has been questioned in different ways. Piddington |H9|, EDI pointed out that small- 



scale magnetic fields can grow large enough (until equipartition is reached) to swamp 
the dynamo action. The quenching of the dynamo action has been numerically shown 
by Kulsrud and Anderson More recently, it has been argued that if the large- 

scale magnetic field reaches the critical value Re^^^v f\ the dynamo action could also be 
quenched [|^, |5^ . 



MHD equations can be derived from a microscopic (kinetic) approach and also from 
a macroscopic approach where the displacement current is neglected |Q . If the displace- 
ment current is neglected the electric field can be expressed using the Ohm law and the 
magnetic diffusivity equation is obtained 

^ = V X (V X 5) + -V'^. (3.6) 
or] a 



The first term at the r.h.s. of Eq. ( |3.6| ) is related to the dynamo term, while the second 
term is the diffusivity term. The conductivity a appearing in Eq. (|3.6| ) is a global quantity 
which can be computed in a kinetic approach during a given phase of evolution of 



the background geometry [Q]. If the plasma is non-relativistic a oc 

If plasma 

is relativistic a oc T. In Eq. ( |3.6| ) the contribution containing the conductivity is usually 
called magnetic diffusivity term. The magnetic diffusivity scale is defined as 

= \[^, (3.7) 
V a 

where tu is the age of the Universe at the corresponding epoch. Typical galactic values 
are 1/a ~ lO^^cmVsec, Hq ~ r^^ ~ IQ-^^ Hz, ~ A.U. . 

Eq. ( p.6|) is exact, in the sense that both V and B contain long and short wavelength 
modes. The aim of the various attempts of the dynamo theory is to get an equation 
describing only the "mean value" of the magnetic field. To this end the fisrt step is to 
separate the exact magnetic and velocity fields as 

B={B) + b, 

V={V)+ V, (3.8) 

where (B) and {V) are the averages over an ensemble of many realizations of the velocity 
field V. In order to derive the standard form of the dynamo equations few important 
assumptions should be made. These assumptions are : 



^ ReM is the magnetic Reynolds number i.e. , approximately, the ratio of the first (the dynamo 
contribution) over the second term (the magnetic diffusivity term) appearing at the right hand side of 



Eq. (3.6); v is the velocity field at the outer scale of turbulence. 
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- The scale of variation of the turbulent motion v should be smaller than the typical 
scale of variation of {B). In the galactic problem (V^) is the differential rotation 
of the galaxy, while v is the turbulent motion generated by stars and supernovae. 
Typically the scale of variation of v is less than 100 pc while the interesting scales 
for {B) are larger than the kpc. 

- The field 6 is such that \h\ < 

- It should happen that ({/ ■ V x {7) 7^ 0. 

- Magnetic flux is frozen into the plasma (i.e. magnetic flux is conserved). 



From the magnetic diffusivity equation ( p.6|) , and using the hsted assumptions, it is pos- 
sible to derive the typical structure of the dynamo term by carefully averaging over the 
velocity field according to the procedure outlined in |l], ||, ^ . Inserting Eq. ( p.8|) into (|3.6|) 
and breaking the equation into a mean part and a random part, two separate induction 
equations can be obtained for the mean and random parts of the magnetic field 



d{B) 
dr] 



V X (^(V) X (5)) + V X (tJx 6), (3.9) 

db ^ ^ ^ ^ ^ ^ ^ 

= V X (t; X (B)) + V X ((V) X b) + V X (v X b) - V X (v X b), (3.10) 

Of] 



where the (magnetic) diffusivity terms have been neglected. In Eq. ( p.9|) , {v x b) is called 
"turbulent emf " and it is the average of the cross product of the small-scale velocity field 
V and of the small scale magnetic field b over a scale much smaller than the scale of {B) 
but much larger than the scale of turbulence. Sometimes, the calculation of the effect 
of {v X b) is done in the case of incompressible and isotropic turbulence. In this case 
{v X b) = 0. This estimate is, however, not realistic since (B) is not isotropic. More 
correctly |3^ , {v x b) should be evaluated by using Eq. (|3.10| ) which is usually written in 



a simplified form 

^ = Vx{vx {B)), (3.11) 

where all but the first term of Eq. ( |3.10| ) have been neglected. To neglect the term 
V X ((V^) X b) does not pose any problem since it corresponds to choose a reference frame 
where {V) is constant. However, the other terms, neglected in Eq. (|3.11| ), are dropped 



because it is assumed that \b\ ^ \{B)\. This assumption may not be valid all the time 
and for all the scales. The validity of Eq. ( p.ll|) seems to require that l/a is very large 



so that magnetic diffusivity can keep always b small |Q . On the other hand |^ one can 
argue that b is only present over very small scales (smaller than 100 pc) and in this case 
the approximate form of eq. ( |3.11|) seems to be more justified. 
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From Eqs. ( p.9|) -( pni) it is possible to get to the final result for the evolution equation 
of {B) 1^ as it is usally quoted 



drj 

where 



V X + (3V'{B) + V X (V) X (fi) , (3.12) 



a = -^{v-V X v), (3.13) 

/3=j{v'), (3.14) 

where a is the dynamo term, f3 is the diffusion term and r is the typical correlation time 
of the velocity field. The term a is, in general, space-dependent. 

It is interesting to point out that the dynamo term in Eq. ( p.l2|) has a simple 
electrodynamical meaning (when a is constant), namely, it can be interpreted as a mean 
ohmic current directed along the magnetic field: 

J = -aB. (3.15) 

This equation tells us that an ensemble of screw-like vortices with zero mean helicity is 
able to generate loops in the magnetic flux tubes in a plane orthogonal to the one of the 
original field. 

If the velocity field is parity-invariant (i.e. no vorticity for scales comparable with 
the correlation length of the magnetic field), then the dynamics of the infrared modes is 
decoupled from the velocity field since, over those scales, a = 0. When the (averaged) 
dynamo term dominates in Eq. ( p.l2|) , magnetic fields can be exponentially amplified. 
The standard lore is that the dynamo action stops when the value of the magnetic field 
reaches the equipartition value (i.e. when the magnetic and kinetic energy of the plasma 
are comparable). At this point the dynamo "saturates". This statement means that, 
in more dynamical terms |l|], back-reaction effects cannot be neglected anymore and Eq. 
( ^.12| ) should then be supplemented with non-linear terms (of the order of B"^), whose 
effect is to stabilize the amplification of the magnetic field. 

It is sometimes useful to recall that the full MHD equations can be studied in two 
different limits : the ideal (or superconducting) approximation where the conductivity 
is assumed to be very high and the real (or resistive) limit where the conductivity takes 
a finite value. In the ideal limit both the magnetic flux and the magnetic helicity are 



conserved. This means, formally |55| 



/ B-d^ = -- I V X V X 5 ■ dE, (3.16) 



dr] Jt. a Jt. 

where S is an arbitrary closed surface which moves with the plasma. If we are in the in 
the inertial regime (i.e. L > L^r) we can say that the expression appearing at the right 
hand side is sub- leading and the magnetic flux lines evolve glued to the plasma element. 
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The other quantity which is conserved in the superconducting hmit is the magnetic 
hehcity 

Hm = I d^xA ■ B, (3.17) 



V 



where A is the vector potential []. In Eq. (|3.17|) the vector potential appears and, therefore 



it might seem that the expression is not gauge invariant. This is not the case. In fact 
A- B is not gauge invariant but, none the less, Hm is gauge- invariant since the integration 
volume is defined in such a way that the magnetic field B is parallel to the surface which 
bounds V and which we will call dV. If n is the unit vector normal to dV then B ■ n = 
on dV [|121 . 

The magnetic gyrotropy 

B-VxB (3.18) 
it is a gauge invariant measure of the diffusion rate of Hm at finite conductivity. In fact 

-^Hm = -- I d^x B-V X B. (3.19) 
dr] o Jv 

The magnetic gyrotropy is a useful quantity in order to distinguish different mechanisms 
for the magnetic field generation. Some mechanisms are able to produce magnetic fields 
whose flux lines have a topologically non-trivial structure (i.e. (i? ■ V x 5) 7^ 0). 

Usually the picture for the formation of galactic magnetic fields is related to the 
possibility of implementing the dynamo mechanism. By comparing the rotation period 
with the age of the galaxy (for a Universe with ~ 0.7, h ~ 0.65 and ~ 0.3) 
the number of rotations performed by the galaxy since its origin is approximately 30. 
During these 30 rotations the dynamo term of Eq. ( p.l2|) dominates against the magnetic 
diffusivity. As a consequence an instability develops. This instability can be used in order 
to drive the magnetic field from some small initial condition up to its observed value. 
Eq. ( p.l2| ) is linear in the mean magnetic field. Hence, initial conditions for the mean 
magnetic field should be postulated at a given time and over a given scale. This initial 



mean field, postulated as initial condition of (|3.12|) is usually called seed. 

Most of the work in the context of the dynamo theory focuses on reproducing the 
correct features of the magnetic field of our galaxy. The achievable amplification produced 
by the dynamo instability can be at most of 10^^, i.e. e^°. Thus, if the present value of 
the galactic magnetic field is 10~^ Gauss, its value right after the gravitational collapse of 
the protogalaxy might have been as small as 10~^^ Gauss over a typical scale of 30-100 
kpc. 

There is a simple way to relate the value of the magnetic fields right after gravitational 
collapse to the value of the magnetic field right before gravitational collapse. Since the 
gravitational collapse occurs at high conductivity the magnetic flux and the magnetic 



^Notice that in conformally flat FRW spaces the radiation gauge is conformany invariant. This 
property is not shared by the Lorentz gauge condition |5q|. 
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helicity are both conserved. Right before the formation of the galaxy a patch of matter 
of roughly 1 Mpc collapses by gravitational instability. Right before the collapse the 
mean energy density of the patch, stored in matter, is of the order of the critical density 
of the Universe. Right after collapse the mean matter density of the protogalaxy is, 
approximately, six orders of magnitude larger than the critical density. 

Since the physical size of the patch decreases from 1 Mpc to 30 kpc the magnetic field 
increases, because of flux conservation, of a factor (pa/Pb)^''^ ~ 10"^ where pa and pb are, 
respectively the energy densities right after and right before gravitational collapse. The 
correct initial condition in order to turn on the dynamo instability is 5 ~ 10~^^ Gauss 
over a scale of 1 Mpc, right before gravitational collapse. 

Since the flux is conserved the ratio between the magnetic energy density, Pb{L, rf) and 
the energy density sitting in radiation, p.^{ri) is almost constant and therefore, in terms 
of this quantity (which is only scale dependent but not time dependent), the dynamo 
requirement can be rephrased as 

r^{^L) = P^^^>lQ-^\ L~lMpc, (3.20) 

to be compared with the value tb ~ 10^^ which would lead to the galactic magnetic field 
only thanks to the collapse and without the need of dynamo action. This is the case when 
the magnetic field is fully primordial. 

The estimate of Eq. ( |3.20| ) is, to say the least, rather generous and has been presented 
just in order to make contact with several papers (concerned with the origin of large scale 
magnetic fields) using such an estimate. Eq. ( |3.20D is based on the (highly questionable) 
assumption that the amplification occurs over thirty e-folds while the magnetic flux is 
completely frozen in. In the real situation, the achievable amplification is much smaller. 
Typically a good seed would not be 10~^^ G after collapse (as we assumed for the simplicity 
of the discussion) but of the order of (or larger than) 10~^^-10~^^ G . 



The possible applications of dynamo mechanism to clusters is still under debate and 
it seems more problematic [|l^, ^ The typical scale of the gravitational collapse 
of a cluster is larger (roughly by one order of magnitude) than the scale of gravitational 
collapse of the protogalaxy. Furthermore, the mean mass density within the Abell radius ( 
~ 1.5/1"^ Mpc) is roughly 10^ larger than the critical density. Consequently, clusters rotate 
less than galaxies since their origin and the value of t-q^L) has to be larger than in the case 
of galaxies. Since the details of the dynamo mechanism applied to clusters are not clear, 
at present, it will be required that rB(i^Mpc) ^ 10~^^ [for instance rB(i^Mpc) — 10^^^]. 



4 A twofold path to the origin 

Back in the late sixties harrison [^] suggested that the initial conditions of the magnetic 
diffusivity equation might have something to do with cosmology in the same way as 
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he suggested that the primordial spectrum of gravitational potential fluctuations (i.e. 
the Harrison- Zeldovich spectrum) might be produced in some primordial phase of the 
evolution of the Universe. Since then, several mechanisms have been invoked in order to 
explain the origin of the magnetic seeds |62, Qi\ and few of them are compatible with 



inflationary evolution. It is not my purpose to review here all the different mechanisms 
which have been proposed and good reviews exist already |61|. A very incomplete 



selection of references is, however, reported ||6^, |6^, Furthermore, more details on 
this topic can be found in the contribution of D. Boyanovsky |l65| . 

In spite of the richness of the theoretical models, the mechanisms for magnetic field 
generation can be divided, broadly speaking, into two categories: astrophysical ^ 



and cosmological. The cosmological mechanisms can be divided, in their turn, into causal 
mechanisms (where the magnetic seeds are produced at a given time inside the horizon) 
and inflationary mechanisms where correlations in the magnetic field are produced outside 
the horizon. Astrophysical mechanisms have always to explain the initial conditions of 
Eq. ( |3.12| ). This is because the MHD are linear in the magnetic fields. It is questionable 
if purely astrophysical considerations can set a natural initial condition for the dynamo 
amplification. 



4.1 Turbulence? 



Causal mechanisms usually fail in reproducing the correct correlation scale of the field 
whereas inflationary mechanisms have problems in reproducing the correct amplitude 
required in order to turn on successfully the dynamo action. In the context of causal 
mechanisms there are interesting proposals in order to enlarge the correlation scale. These 
proposals have to do with the possible occurrence of turbulence in the early Universe. The 
ratio of the magnetic Reynolds number to the kinetic Reynolds number is the Prandtl 
number 

Rgm 



Pr 



M 



Re 



(4.1) 



where u is the thermal diffusivity coefficient and a, is, as usual, the conductivity. Consider, 
for instance, the case of the electroweak epoch ^ 0. At this epoch taking 

~ 3cm we get that Ptm ua ^ 10^ where the bulk velocity of the plasma is of the 
order of the bubble wall velocity at the epoch of the phase transition. 

This means that the early universe is both kinetically and magnetically turbulent. The 
features of magnetic and kinetic turbulence are different. This aspect reflects in a spectrum 
of fluctuations is different from the usual Kolmogorov spectrum |^|. If the Universe is 
both magnetically and kinetically turbulent it has been speculated that an inverse cascade 
mechanism can occur [BB, B^, 17^. This idea was originally put forward in the context 



of MHD simulations The inverse cascade would imply a growth in the correlation 
scale of the magnetic inhomogeneities and it has been shown to occur numerically in 



the approximation of unitary Prandtl number [43|. Specific cascade models have been 
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also studied [g^, ^ A particularly important role is played, in this context, by 
the initial spectrum of magnetic fields (the so-called injection spectrum |6^) and by the 
topological properties of the magnetic fiux lines. If the system has non vanishing magnetic 
helicity and magnetic gyrotropy it was suggested that the inverse cascade can occur more 
efficiently ||73|, Recently simulations have discussed the possibility of inverse cascade 
in realistic MHD models |]73|. More analytic discussions based on renormalization group 
approach applied to turbulent MHD seem to be not totally consistent with the occurrence 



of inverse cascade at large scales 75 



4.2 Magnetic fields from dynamical gauge couplings 

Large scale magnetic fiuctuations can be generated during the early history of the Universe 
and can go outside the horizon with a mechanism similar to the one required in order to 
produce fiuctuations in the gravitational (Bardeen) potential. In this case the correlation 
scale of the magnetic inhomogeneities can be large. However, the typical amplitudes 
obtainable in this class of models may be too small. 

The key property allowing the amplification of the fiuctuations of the scalar and tensor 
modes of the geometry is the fact that the corresponding equations of motion are not 
invariant under Weyl rescaling of a (conformally fiat) metric of FRW type. In this sense 
the evolution equations of relic gravitons and of the scalar modes of the geometry are said 
to be not conformally invariant. If gauge couplings are dynamical, the evolution equations 
of the gauge field are also not conformally invariant. Interesting examples in this direction 
are models containing extra-dimensions and scalar-tensor theories of gravity where the 
gauge coupling is, effectively, a scalar degree of freedom evolving in a given geometry. 

The remarkable similarity of the abundances of light elements in different galaxies 
leads to postulate that the Universe had to be dominated by radiation at the moment 
when the light elements were formed, namely for temperatures of approximately 0.1 MeV 
||76| , [77| . Prior to the moment of nucleosynthesis even indirect informations concerning 
the thermodynamical state of our Universe are lacking even if our knowledge of particle 
physics could give us important hints concerning the dynamics of the electroweak phase 



transition 7S 



The success of big-bang nucleosynthesis (BBN) sets limits on alternative cosmological 
scenarios. Departures from homogeneity and isotropy of the background geometry 



can be successfully constrained. In the same spirit, BBN can also set limits on the 
dynamical evolution of internal dimensions |5^. Internal dimensions are an essential 



ingredient of theories attempting the unification of gravitational and gauge interactions in 
a higher dimensional background like Kaluza-Klein theories and superstring theories 



84| 



Defining, respectively, Iibbn and ho as the size of the internal dimensions at the BBN 
time and at the present epoch, the maximal variation allowed to the internal scale factor 



from the BBN time can be expressed as hBBN /bo ^ 1 + e where |e| < 10 ||8T|, |8^. The 
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bounds on the variation of the internal dimensions during the matter dominated epoch 
are even stronger. Denoting with an over-dot the derivation with respect to the cosmic 
time coordinate, we have that \b/h\ < lO^^ifo where Hq is the present value of the Hubble 



parameter The fact that the time evolution of internal dimensions is so tightly 

constrained for temperatures lower of 1 MeV does not forbid that they could have been 
dynamical prior to that epoch. Moreover, recent observational evidence [^, |8^, |87 



seem 



to imply that the fine structure constant can be changing even today. 

Suppose that prior to BBN internal dimensions were evolving in time and assume, 
for sake of simplicity, that after BBN the internal dimensions have been frozen to their 
present (constant) value. Consider a homogeneous and anisotropic manifold whose line 
element can be written as 

ds^ = G^ydx^dx'^ = a?{7])[drf' — '^ijdx^dx^] — b'^{ri)^abdy"'dy^, 
^,1/ = 0, D — 1 = d + n, i,j = l,...,d, a,b = d + 1, ...,d + n. (4.2) 

[rj is the conformal time coordinate related, as usual to the cosmic time t = J a{ri)dri 
; 7jj(x), 'Jabiy) are the metric tensors of two maximally symmetric Euclidean manifolds 
parameterized, respectively, by the "internal" and the "external" coordinates {x*} and 
{y'^}]. The metric of Eq. ([4.2|) describes the situation in which the d external dimensions 
(evolving with scale factor a{T])) and the n internal ones (evolving with scale factor b{ri)) 



are dynamically decoupled from each other The results of the present investigation 



however, can be easily generalized to the case of n different scale factors in the internal 
manifold. 

Consider now a pure electromagnetic fluctuation decoupled from the sources, repre- 
senting an electromagnetic wave propagating in the d-dimensional external space such 
that Afj, = Af^{x,ri), Aa = 0. In the metric given in Eq. (|4.2|) the evolution equation of 
the gauge field fluctuations can be written as 

^ ^J^/^G''^G'"'F^^=0, (4.3) 



where Fa/3 = V^qA^j is the gauge field strength and G is the determinant of the D di- 
mensional metric. Notice that if n = the space-time is isotropic and, therefore, the 
Maxwell's equations can be reduced (by trivial rescaling) to the flat space equations. If 
n 7^ we have that the evolution equation of the electromagnetic fluctuations propagat- 
ing in the external rf-dimensional manifold will receive a contribution from the internal 
dimensions which cannot be rescaled away. 

In the radiation gauge {Aq = and VjA* = 0) the evolution the vector potentials can 
be written as ^ 

A'^ + nTJ^, - V'A, = 0, ^ =\- (4.4) 

The vector potentials Ai are already rescaled with respect to the (conformally flat) (i + 1 
dimensional metric. In terms of the canonical normal modes of oscillations Ai = b'^^'^Aj 
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the previous equation can be written in a simpler form, namely 



- V{r])A^ - VMi = 0, 



(4.5) 



From this set of equations the induced large scale magnetic fields can be computed in var- 



ious models for the evolution of the internal manifold p9[ . It should be noticed that large 
magnetic seeds are produced in this context only if internal dimensions are rather large 
if compared to the Planck length. This requires a careful discussion of the localization 
properties of gauge fields in the presence of large extra-dimensions |^ which is beyond 
the scope of the present discussion. 

The same effect of magnetic field generation can be also obtained in the case of time- 
evolving gauge coupling already in four dimensions. In order to emphasize this phe- 
nomenon it will now be shown how squeezed states of relic photons can be produced 
9T|. It will be imagined that quantum- mechanical fluctuations of the gauge field will be 



present at some initial stage in the evolution of the Universe. 

The squeezed states formalism has been successfully applied to the analysis of tensor, 
scalar |92| and rotational fluctuations of the metric by Grishchuk and collaborators. In 
the case of relic gravitons and relic phonons the analogy with quantum optics is certainly 
very inspiring. In the case of relic photons the analogy is even closer since the time 
variation of the dilaton coupling plays directly the role of the laser "pump" which is 



employed in order to produce experimentally observable squeezed states [94 



The effective action of a generic Abelian gauge field in four space-time dimensions 
reads 



(4.6) 



where F^fs = V[q,A^] is the Maxwell field strength and Vq, is the covariant derivative with 
respect to the string frame metric G^u- In Eq- ( |4.6|) g is the (four dimensional) gauge 
coupling which is related to the expectation value of a scalar degree of freedom. 

From Eq. ( [4. 61) it is possible to derive the Hamiltonian and the Hamiltonian density 
of the gauge field fluctuations 



+e{g)d_k,adk,a + e*(5')4,aO-fc 



k (dl ^dk,a + alfc ad-k,a + 1) 



id) 
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(4.7) 



where The (two-modes) Hamiltonian contains a free part and the effect of the variation of 
the coupling constant is encoded in the (Hermitian) interaction term which is quadratic in 
the creation and annihilation operators whose evolution equations, read, in the Heisenberg 
picture 



ddk,c 
drj 
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The general solution of the previous system of equations can be written in terms of a 
Bogoliubov-Valatin transformation 

^iaiv) = fJ-lMbla + KMb-k,a (4.9) 

where ak,a{0) = bk,a and d-k,a{0) = b_k,a- Unitarity requires that the two complex 
functions fik{v) and I'kiv) are subjected to the condition 1/1^(7]) p — = 1 which also 

implies that fik{v) and i^kiv) can be parameterized in terms of one real amplitude and 
two real phases 

fj^^ = e^^fe cosh Tfc, Uk = 6*^^'^'="^'=) sinh r^, (4.10) 

(r is sometimes called squeezing parameter and (pk is the squeezing phase; from now on 
we will drop the subscript labeling each polarization if not strictly necessary). The total 
number of produced photons 

{0-kOk\al{r])ak{r]) + dl^.a-k\OkO-k) = 2 Uk. (4.11) 

is expressed in terms of = sinh^r^, i.e. the mean number of produced photon pairs in 
the mode k. Inserting Eqs. ([4.9|),( ^rT0D and ( [4.11| ) into Eqs. ( [4.8| ) we can derive a closed 
system involving only the Uk and the related phases: 

^ = -2/(n,,)-cos20fc, (4.12) 

— = -A; + — — — sm20fc, (4.13) 
dv g f{nk) 

"^^'^ , , g'dfijh) . ^ , , .^ 

—k-\ — — sm 2(pk, (4.14) 



drj g duk 

where f(nk) = \Jnk{nk + 1). 

The two-point function of the magnetic fields 

Gij{r,Tl) = (0_feOfc|A(^,^)4(^+^^)|OfeO-fc) (4.15) 
can be expressed, using Eqs. ( [4.9D and (|4.1CI|) 

e,(.-)=/Ac;«We-- (4.16) 

where 

Qij{k,7]) = ^^^--y^-j^A;[2sinhVfc + sinh 2rfc cos 20fc] 

= Y.<{k)e-{k) = (5,, - (4.17) 
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(the vacuum contribution, occurring for = 0, has been consistently subtracted). The 
intercept for f= of the two-point function traced with respect to the two polarizations 
is related to the magnetic energy density 
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■ [2 sinh^ Tfc + sinh 2rfc cos 20fc] (4.18) 



(where a; = A;/a is the physical frequency). The two-point function and its trace only 
depend upon Uk and upon 0^. Since Eqs. ( [4.12| ) and ( [4.14[ ) do not contain any dependence 



upon 9k we can attempt to solve the time evolution by solving them simultaneously. In 
terms of the new variable x = krj Eqs. ( [4. 121 ) and ( |4.14| ) can be written as 



— = -l + — -^sm20fc, (4.19) 

ax ax drik 

-2/(rrfc)cos20fc, (4.20) 



ding 

If \{d\ng/dx){df(nk)/dnk) sin 20^1 > 1, then Eqs. (|4.19|) and (|4.20| ) can be written as 



du, _ ^dj^^^^^ ^ -2/(,..)i^ (4.21) 



d\ng drik d\n.g 1 + ui. 

where (pk = arctanwfc. By trivial algebra we can get a differential relation between Uk 
and Uk which can be exactly integrated with the result that m| — f(nk)uk + 1 = 0. 
By inverting this last relation we obtain two different solutions with equivalent physical 
properties, namely 

Mnk) = [^{\Jnk(nk + 1) ± ^nkijik + 1) - 4)]. (4.22) 

If we choose the minus sign in Eq. ( |4.22| ) we obtain that 0^ ~ (m + l)vr/2, m = 0, 1, 2... 
with corrections of order l/rik . In the opposite case (pk ~ arctan (rIfc/2) within the 
same accuracy of the previous case(i.e. l/n^)- By using the relation between Uk and 
rik the condition \{d\ng/dx){df(jik)/dnk) sin 2(j)k\ > 1 is equivalent to x^l, if, as we 
are assuming, \g'/g\ vanishes as ?7~^ for t] ±oo and it is, piece-wise, continuous. 
By inserting Eq. (|4.22| ) into Eq. ([4.19|) a consistent solution can be obtained, in this 



case, if we integrate the system between r]f and r]i defined as the conformal times where 
\{dhig/dx){df(nk)/dnk)sm2(l)k\ = 1: 



- ( \ ^ (aim) gjvi) \ 

"'<"■''- 31^ -S(^;Tj 



TT 1 

2+ ^^)- 



()fc(r/) ~ (m + 1)^ + 0(^1^), m = 0,l,2... (4.23) 
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If \{d\n. g / dx){df {rik) / duk) sin 20^1 < 1 (i. e. a; > 1) the consistent solution of our system 
is given by 



nk{r]f) = sinh^ ^2 j \YLg{x')sm.2x' dx' 
(pk ~ —krj + v^fc, — constant. 



(4.24) 



If the coupling constant evolves continuously between — oo and +oo with a (global) maxi- 
mum located at some time rj^. then, for a; > 1, 71^ ~ const.. Indeed by taking trial functions 
with bell-like shape for \g' / g\ we can show that oscillates around zero for large (pk- 

From Eq. ( [4.23| ) the magnetic energy density of Eq. ([4.18|) can be computed in 
different scenarios and related to the ratio discussed in Eq. ( p.20[ ). Recalling that the 
present frequency corresponding to 1 Mpc is roughly uoq ~ 10~^^ Hz, the ratio r{ujQ) can 
be estimated. Time evolution of gauge coupling during an inflationary phase produces 
rather large seeds r{ujG) ~ 10~^^ |9^. Furthermore, pre-big bang models also lead 
to large seeds which can be even r{ujQ) ~ 10~* |100| |. There is a difference in the 
models discussed in and in |10(]|| . In the case of the gauge coupling is related 
to the expectation value of a scalar field which evolves during an inflationary phase of 
de Sitter (or quasi de Sitter type). This scalar degree of freedom is not the inflaton and 



it is not a source of the background geometry. On the contrary, in ||99|, |100|| the gauge 
coupling is a source of the evolution of the background geometry since it is connected to 
the expectation value of the dilaton field whose specific evolution dictates the nature of 
the pre-big bang solutions used in order to describe the dynamics of the Universe in its 
early stages. 



5 EWPT and BAU 

In the previous Section it has been stressed that causal mechanisms have, in general 
problems with the correlation scale of the obtained field, while inflationary mechanisms 
may have problems with the seed amplitude. In spite of this, it should be borne in mind 
that magnetic fields are generated over all physical scales compatible with the plasma 
dynamics at a given epoch. Hence, even if the magnetic fields at large scales may be 
very minute, magnetic fields at smaller scales may have a very interesting impact on 
different moments of the life of the Universe. The physical picture we have in mind is 
then the following. Suppose that conformal invariance is broken at some stage in the 
evolution of the Universe, for instance thanks to the (effective) time variation of gauge 
couplings. Then, vacuum fluctuations will go outside the horizon and will be amplified. 
The amplified magnetic inhomogeneities will re-enter (crossing the horizon a second time) 
during different moments of the life of the Universe and, in particular, even before the 
BBN epoch. 

In the following various effects of these magnetic fields will be considered starting with 
the EW epoch. The electroweak epoch occurs when the temperature of the plasma was 
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roughly T Tc ^ 100 GeV. The physical size of the horizon was, at that time , H^^ ~ 3 
cm. The electroweak epoch occurs, approximately, when the Universe was 10~^^ sec old. 

At small temperatures and small densities of different fermionic charges the SUl{'2) ® 
?7y(l) is broken down to the f/em(l) and the long range fields which can survive in the 
plasma are the ordinary magnetic fields. However, for sufficiently high temperatures the 
SUl{2) ® Uy{1) is restored and non-screened vector modes correspond to hypermagnetic 
fields. At the electroweak epoch the typical size of the horizon is of the order of 3 cm . 
The typical diffusion scale is of the order of 10~^ cm. Therefore, over roughly eight orders 
of magnitude hypermagnetic fields can be present in the plasma without being dissipated 
|70| . The evolution of hypermagnetic fields can be obtained from the anomalous magne- 



tohydrodynamical (AMHD) equations. The AMHD equations generalize the treatment 
of plasma effects involving hypermagnetic fields to the case of finite fermionic density ||7T 



Depending on their topology, hypermagnetic fields can have various consequences 



ffOl , [71|| . If the hypermagnetic flux lines have a trivial topology they can have an im- 
pact on the phase diagram of the electroweak phase transition | |1 1 1| , |112|| . If the topology 
of hypermagnetic fields is non trivial, hypermagnetic knots can be formed ||104|| and, un- 



der specific conditions, the BAU can be generated | |105| |. The gauge field fluctuations 



produced as a result of the parametric amplification of vacuum fluctuations always lead 
to hypermagnetic fields with topologically trivial structure (i.e. with zero magnetic he- 
licity and gyrotropy). However, thanks to pseudo-scalar couplings, a topologically trivial 
background of hypermagnetic flux lines may lead to a non-zero magnetic gyrotropy and, 
hence, to some kind of hypermagnetic knot with topologically non-trivial structure. 

A classical hypermagnetic background in the symmetric phase of the EW theory can 
produce interesting amounts of gravitational radiation in a frequency range between 10"'^ 
Hz and the kHz. The lower tail falls into the LISA window while the higher tail falls in 
the VIRGO/LIGO window. For the hypermagnetic background required in order to seed 
the BAU the amplitude of the obtained GW can be even six orders of magnitude larger 
than the inflationary predictions. In this context, the mechanism of baryon asymmetry 
generation is connected with GW production ||104|, |105|. 



5.1 Hypermagnetic knots 



It is possible to construct hypermagnetic knot configurations with finite energy and he- 
licity which are localized in space and within typical distance scale Lg. Let us consider in 
fact the following configuration in spherical coordinates [105| 



2Bn COS 9 



TtLs [7^2 + 1] 

2Bo sine 



2 ' 
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(5.1) 



where TZ = r / Lg is the rescaled radius and Bq is some dimensionless amphtude and n is 
just an integer number whose physical interpretation will become clear in a moment. The 
hypermagnetic field can be easily computed from the previous expression and it is 



ABq n cos 9 
45o 7^2 - 1 
8Bo 7^sin^ 



r?T,sin6', 



(5.2) 



The poloidal and toroidal components of H can be usefully expressed as Hp = Hrer+Heee 
and Ht = H-^e^. The Chern-Simons number is finite and it is given by 
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/2 



CS 



y ■ Hyd^X 



g'2 roo ^nBl 7^2rf7^ g'^nBl 



3271^ Jv^ '"-^ 327r2 7o vr^ [7^2 + 1]^ 32n^ 

We can also compute the total helicity of the configuration namely 



(5.3) 



- - - ,o 256 52 n U^dU 

Hy-Vx nyd'x = -f- / 

V ttL^ jo {l + 7i^)^ 



5B^n 



(5.4) 



We can compute also the total energy of the field 

1 



E 



2Jv ' ' 2 L 



(5.5) 



and we discover that it is proportional to n^. This means that one way of increasing the 
total energy of the field is to increase the number of knots and twists in the fiux lines. We 
can also have some real space pictures of the core of the knot (i.e. 71 = r/Lg < 1). This 
type of configurations can be also obtained by projecting a non-Abelian SU(2) (vacuum) 
gauge field on a fixed electromagnetic direction ||106|| f\ These configurations have been also 
studied in [[108| , p.09|| . In particular, in [p.09|] , the relaxation of HK has been investigated 
with a technique different from the one employed in ||104| , |105|| but with similar results. 
The problem of scattering of fermions in the background of hypermagnetic fields has been 



also studied in 



II In order to interpret these solutions it is very interesting to make use of the Clebsh decomposition. 
The imphcations of this decomposition (beyond the hydrodynamical context, where it was originaUy dis- 
covered) have been recently discussed (see [lO?] and references therein). I thank R. Jackiw for interesting 
discussions about this point. 
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Topologically non-trivial configurations of tlie liypermagnetic flux lines lead to the 
formation of liypermagnetic knots (HK) whose decay might seed the Baryon Asymmetry 
of the Universe (BAU) . HK can be dynamically generated provided a topologically trivial 
(i.e. stochastic) distribution of flux lines is already present in the symmetric phase of 
the electroweak (EW) theory ||104| , |105|| . In spite of the mechanism generating the HK, 



their typical size must exceed the diffusivity length scale. In the minimal standard model 
(MSM) (but not necessarily in its supersymmetric extension) HK are washed out. 

The importance of the topological properties of long range (Abelian) hj^ercharge 



magnetic fields has been stressed in the past [p.l3| , p.l4| , p.l5| , p.l6|| . In [ p.l7|| it was argued 
that if the spectrum of hypermagnetic fields is dominated by parity non-invariant Chern- 
Simons (CS) condensates, the BAU could be the result of their decay. Most of the 
mechanisms often invoked for the origin of large scale magnetic fields in the early Universe 
seem to imply the production of topologically trivial (i.e. stochastic) configurations of 



magnetic fields |62 



5.2 Hypermagnetic knots and BAU 

Suppose that the EW plasma is filled, for T > with topologically trivial hypermagnetic 
fields Tiy, which can be physically pictured as a collection of fiux tubes (closed because of 
the transversality of the field lines) evolving independently without breaking or intersect- 
ing with each other. If the field distribution is topologically trivial (i.e. (?iyVx?-^y) = 0) 
parity is a good symmetry of the plasma and the field can be completely homogeneous. 
We name hypermagnetic knots those CS condensates carrying a non vanishing (averaged) 
hypermagnetic helicity (i.e. (TYy ■ V x Tiy) 7^ 0). If (?iy ■ V x Tiy) 7^ parity is broken 
for scales comparable with the size of the HK, the fiux lines are knotted and the field Tiy 
cannot be completely homogeneous. 

In order to seed the BAU a network of HK should be present at high temperatures 



||70| , [TTi |117|| . In fact for temperatures larger than the fermionic number is stored both 
in HK and in real fermions. For T < Tc, the HK should release real fermions since the 
ordinary magnetic fields (present after EW symmetry breaking) do not carry fermionic 
number. If the EWPT is strongly first order the decay of the HK can offer some seeds for 
the BAU generation ||117|| . This last condition can be met in the minimal supersymmetric 



standard model (MSSM) 1TTB|, [11^1, pq, pi 



Under these hypotheses the integration of the f/(l)y anomaly equation ||117|| gives the 
CS number density carried by the HK which is in turn related to the density of baryonic 
number Ub for the case of fermionic generations. 

ns, , a' nf{nY-^xnY)MoT , g'^ , ^, 

— [tc) = 7, TT-^ ^ =ir 

s 27r(Tc s T + T-}i An 

{g' is the f/(l)y coupling and s = (2/45)7r^A''e//T^ is the entropy density; Nejf, at Tc, 
is 106.75 in the MSM; Mq = Mp/lM^/NJf ~ 7.1 x lO^^GeV). In Eq. T is 
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the perturbative rate of the right electron chirahty flip processes (i.e. scattering of right 
electrons with the Higgs and gauge bosons and with the top quarks because of their large 
Yukawa coupling) which are the slowest reactions in the plasma and 

783 \ny\' 

is the rate of right electron dilution induced by the presence of a hypermagnetic field. In 
the MSM we have that F < Tj-c ||122|| whereas in the MSSM F can naturally be larger than 
Ft.^. Unfortunately, in the MSM a hypermagnetic field can modify the phase diagram of 
the phase transition but cannot make the phase transition strongly first order for large 



masses of the Higgs boson ||111|| . Therefore, we will concentrate on the case F > T-^ 
and we will show that in the opposite limit the BAU will be anyway small even if some 
(presently unknown) mechanism would make the EWPT strongly first order in the MSM. 



HK can be dynamically generated [|104| , p.05|l (see also [|123|| ). Gauge-invariance and 



transversality of the magnetic fields suggest that perhaps the only way of producing 
(?iy Vx?-^y) 7^ is to postulate, a time-dependent interaction between the two (physical) 
polarizations of the hypercharge field Y^. Having defined the Abelian field strength Yaf3 = 
V[aY/^] and its dual Yap such an interaction can be described, in curved space, by the 
Lagrangian density 



(5. 



where Qf^i, is the metric tensor and g its determinant, c is the coupling constant and M 
is a typical scale. This interaction is plausible if the U{1)y anomaly is coupled, (in the 
symmetric phase of the EW theory ) to dynamical pseudoscalar particles ip (like the axial 
Higgs of the MSSM). Thanks to the presence of pseudoscalar particles, the two polar- 
izations of Tiy evolve in a slightly different way producing, ultimately, inhomogeneous 
HK. 

Suppose that an inflationary phase with a(r) ~ is continuously matched, at the 
transition time Ti, to a radiation dominated phase where a(r) ~ r. Consider then a 
massive pseudoscalar field ip which oscillates during the last stages of the inflationary 
evolution with typical amplitude ipo ~ M. As a result of the inflationary evolution 
\Vip\ -C i^'. Consequently, the phase of ip can get frozen. Provided the pseudoscalar 
mass m is larger than the inflationary curvature scale Hi ~ const., the ip oscillations are 
converted, at the end of the quasi-de Sitter stage, in a net helicity arising as a result of 
the different evolution of the two (circularly polarized) vector potentials 

Y± + (tY^ + ojIY± = 0, Hy = VxY (5.9) 
col = k"^ ^ k—aip (5.10) 
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(where we denoted with Hy = a^Tiy the curved space fields and with a = the rescaled 
hyperconductivity; the prime denotes derivation with respect to conformal time r whereas 
the over-dot denotes differentiation with respect to cosmic time t). 

Since ujj^ ^ uj^ the hehcity gets amphfied according to Eq. (|5.9|) . There are two 
important points to stress in this context. First of all the plasma effects as well as the 
finite density effects are important. This means, in practical terms, that the dissipation 
scales of the problem should be borne in mind. This has not always been done. The second 
point is related to the first. By treating, consistently, the plasma and finite density effects 



in the context of AMHD [^, p. 04] , p.05|| , one realizes that the pseudoscalar coupling of ip 
together with the coupling of the chemical potential are not sufficient in order to seed 
the BAU unless some hypermagnetic background is already present. In other words the 
scenario which leads to the generation of the BAU is the following. Correlation in the 
hypermagnetic fields are generated outside the horizon during infiation by direct breaking 
of conformal invariance. An example in this direction has been given in the previous 
Section. Then hypermagnetic fields re-entering at the electroweak epoch will participate 
in the dynamics and, in particular, they will feel the effect of the anomalous coupling 
either to ip or to the chemical potential. The effect of the anomalous coupling will not be 
to amplify the hypermagnetic background. The anomalous coupling will only make the 
topology of the hypermagnetic flux lines non-trivial. So the statement is that conformal 
invariance is broken and, z/hypermagnetic have anomalous couplings, then a BAU > 10 



-10 



can be achieved without spoiling the standard cosmological evolution [|104|, |105|| . It is 
worth mentioning that this type of scenario may be motivated by the low energy string 
effective action where, by supersymmetry, Kalb-Ramond axions and dilatons are coupled. 



respectively, to Y^^jY^'^ and to the gauge kinetic term ||124| | . It is interesting to notice 



that, in this scenario, the value of the BAU is determined by various particle physics 
parameters but also by the ratio of the hypermagnetic energy density over the energy 
density sitting in radiation during the electroweak epoch, namely, using the language of 
the previous Sections [|104| , |105|| , 

nu , 
— ocr. 5.11 

s 

In order to get a sizable BAU r should be at least 10~^ if the anomalous coupling operates 
during a radiation phase. The value of r could be smaller in models where the anomalous 



coupling is relevant during a low scale inflationary phase ||105 



6 BBN and matter— antimatter fluctuations 

Large scale magnetic flelds possibly present at the BBN epoch can have an impact on the 
light nuclei formation. By reversing the argument, the success of BBN can be used in 
order to bound the magnetic energy density possibly present at the time of formation of 
light nuclei. 
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These bounds are qualitatively different from the ones previously quoted and coming, 
alternatively, from homogeneity [|7^] and isotropy of the background geometry at the 



BBN time. As elaborated in slightly different frameworks through the years ||126| , |12?1 , |12S 



129 , 130], magnetic fields possibly present at the BBN epoch could have a twofold effect. 
On one hand they could enhance the rate of reactions (with an effect proportional to aps) 
and, on the other hand they could artificially increase the expansion rate (with an effect 
proportional to ps)- It turns out that the latter effect is probably the most relevant ||129|| . 
In order to prevent the Universe from expanding too fast at the BBN epoch ps < 0.27pi, 
where p^, is the energy density contributed by the standard three neutrinos for T < 1 
MeV. 

In the previous Section the case of a topologically non-trivial hypermagnetic back- 
ground has been considered (i.e. {Tiy ^Tiv) 7^ 0). In this Section we will instead assume 
that the hypermagnetic background is topologically trivial (i.e. {T-Ly - V ^TLy) = 0). In this 
case, fluctuations in the baryon to entropy ratio will be induced since {{TCy-'V xHy)'^) ^ 0. 
These fluctuations are of isocurvature type and can be related to the spectrum of hyper- 
magnetic fields at the EWPT. Defining as 



A(r,g = J{6(^) (f,tc)5f— ) i^ + r,Q), (6.1) 



the fluctuations in the baryon to entropy ratio at t = tc [j^, the value of A(r, tc) can 
be ralated to the hypermagnetic spectrum which is determined in terms of its amplitude 
^ and its slope e. A physically realistic situation corresponds to the case in which the 
Green's functions of the magnetic hypercharge flelds decay at large distance (i. e. e > 0) 
and this would imply either "blue" ( e > ) or "violet" (e 3> 1) energy spectra. The case of 
"red" spectra (e < 0) will then be left out of our discussion. The flat spectrum corresponds 
to e ^ 1 and may appear quite naturally in string cosmological models |]99| , |100|| . 

Since the fluctuations in the baryon to entropy ratio are not positive deflnite, they will 
induce fluctuations in the baryon to photon ratio, f] at the BBN epoch. The possible effect 
of matter-antimatter fluctuations on BBN depends on the typical scale of of the baryon 
to entropy ratio at the electroweak epoch. Recalling that for T ~ Tc ~ 100 GeV the size 
of the electrowek horizon is approximately 3 cm, fluctuations whose scale is well inside 
the EW horizon at have dissipated by the BBN time through (anti)neutron diffusion. 
The neutron diffusion scale at is 

r^(Tj = 0.3 cm. (6.2) 

The neutron diffusion scale at T = 1 keV is 10^ m, while, today it is 10~^ pc, i.e. of 
the order of the astronomical unit. Matter-antimatter fluctuations smaller than 10^ m 
annihilate before neutrino decoupling and have no effect on BBN. Two possibilities can 
then be envisaged. We could require that the matter-antimatter fluctuations (for scales 
f > ''"n) are small. This will then imply a bound, in the (^,e) plane on the strength of 
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Figure 2: We report the value of the baryon number fluctuations for different parameters 
of the hypermagnetic background ^ and e. 



the hypermagnetic background. In Fig. |^ such an exclusion plot is reported with the 
full line. With the dashed line the bound implied by the increase in the expansion rate 
is reported. Finally with the dot-dashed line the critical density bound is illustrated for 
the same hypermagnetic background. The second possibility is to study the effects of 
large matter-matter domains. These studies led to a slightly different scenario of BBN 
1 1311 , namely BBN with matter-antimatter regions |13|, |13|, |T3|, |13|, |T3|, Q. The 



idea is to discuss BBN in the presence of spherically symmetric regions of anti-matter 
characterized by their radius ta and by the parameter R, i.e. the matter /antimatter ratio. 
Furthermore, in this scenario the net baryon-to-photon ratio, rj, is positive definite and 
non zero. Antimatter domains larger than 10^ m at 1 keV may survive until BBN and 
their dissipation has been analyzed in detail in [|132| , |133| , |134| , |135| , |136| , |137|| . Antimatter 
domains in the range 

10^ m <rA < 10^ m (6.3) 

at 1 keV annihilate before BBN for temperatures between 70 keV and 1 MeV. Since 
the antineutrons annihilate on neutrons, the neutron to proton ratio gets smaller. As a 
consequence, the ^He abundance gets reduced if compared to the standard BBN scenario. 
The maximal scale of matter-antimatter fluctuations is determined by the constraints 
following from possible distortions of the CMB spectrum. The largest scale is of the order 
of 100 pc (today), corresponding to 10^^ m at 1 keV. Suppose that matter-antimatter 
regions are present in the range of Eq. ( |6.3| ). Then the abundance of ^He get reduced. 
The yield of ^He are reported as a function of R, the matter-antimatter ratio and ta. Now, 
we do know that by adding extra-relativistic species the ^He can be increased since the 
Universe expansion gets larger. Then the conclusion is that BBN with matter-antimatter 
domains allows for a larger number of extra-relativistic species if compared to the standard 
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Figure 3: The parameter space of the hypermagnetic background in the case A(r, t^) < 
n^/s for r > (full line). 

BBN scenario. This observation may have implications for the upper bounds on the 
stochastic GW backgrounds of cosmological origin ||138|| since the extra-relativistic species 
present at the BBN epoch can indeed be interpreted as relic gravitons. 



7 GW backgrounds 



If a hypermagnetic background is present for T > Tc, then, as also discussed in [|125|| in the 
context of ordinary MHD, the energy momentum tensor will acquire a small anisotropic 
component which will source the evolution equation of the tensor fluctuations h^^, of the 
metric g^^: 

h!r + 27^/l^ - V^hij = -IGTrGrip. (7.1) 
where r^p is the tensor component of the energy-momentum tensor [|125|| of the hyper- 



magnetic fields. Suppose now, as assumed in that |7i| has constant amplitude and 



that it is also homogeneous. Then as argued in |141| we can easily deduce the critical 



fraction of energy density present today in relic gravitons of EW origin 

^^gw(to) = ^ ^ p,(T,) ^ N^sT^ (7.2) 

Pc 

(zeq = 6000 is the redshift from the time of matter-radiation, equality to the present time 
to). Because of the structure of the AMHD equations, stable hypermagnetic fields will 
be present not only for u^w ~ k^^/a but for all the range cUgw < oj < uj„ where uj„ is the 
diffusivity frequency. Let us assume, for instance, that Tc ~ 100 GeV and A^e// = 106.75. 
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Figure 4: From ||138|| the ^He yield is illustrated in the (i?, ta) plane for rj = 6 x 10 
As the matter /antimatter ratio decreases, we recover the standard "^He yield. 



Then, the (present) values of Wgw is 



(7,3) 



Thus, co'o-(to) ~ lO^cUew Suppose now that Tc ~ 100 GeV; than we will have that ci;ew(^o) 
10^^ Hz. Suppose now, as assumed in ||111||, that 



\n\/T^>0.3. 

This requirement imposes r ~ 0.1-0.001 and, consequently. 



hl^Gw ^ 10"^ - 10 
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(7.4) 



(7.5) 



Notice that this signal would occurr in a (present) frequency range between 10~^ and 
10^ Hz. This signal satisfies the presently available phenomenological bounds on the 
graviton backgrounds of primordial origin. The pulsar timing bound ( which applies for 
present frequencies up ~ 10"*^ Hz and implies h^Qcw < 10~^) is automatically satisfied 
since our hypermagnetic background is defined for 10~^Hz < uj < lO'^Hz. The large 
scale bounds would imply /iq^qw < 7 x 10^^^ but a at much lower frequency (i.e. 10^^^ 
Hz). The signal discussed here is completely absent for frequencies u < cjcw Notice 
that this signal is clearly distinguishable from other stochastic backgrounds occurring at 
much higher frequencies (GHz region) like the ones predicted by quintessential infiation 
142|| . It is equally distinguishable from signals due to pre-big-bang cosmology (mainly 
in the window of ground based interferometers [HB]). The frequency of operation of the 
interferometric devices (VIRGO/LIGO) is located between few Hz and 10 kHz ||143|| . The 
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frequency of operation of LISA is well below the Hz (i.e. 10~^Hz, approximately). In 
this model the signal can be located both in the LISA window and in the VIRGO/LIGO 
window due to the hierarchy between the hypermagnetic diffusivity scale and the horizon 
scale at the phase transition ||104|, |105| . 
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Figure 5: The stochastic background of GW produced by inflationary models with fiat 
logarithmic energy spectrum, illustrated together with the GW background of hypermag- 
netic origin. The frequencies marked with dashed lines correspond to the electroweak 
frequency and to the hypermagnetic diffusivity frequency. 



8 Faraday rotation of CMB ? 



Large scale magnetic fields present at the decoupling epoch can have various consequences. 
For instance they can induce fluctuations in the CMB |^9|, 100], they can distort the 
Planckian spectrum of CMB ||139|| , they can distort the acoustic peaks of CMB anisotropics 



140|| and they can also depolarize CMB ||144 



The polarization of the CMB represents a very interesting observable which has been 
extensively investigated in the past both from the theoretical 
points of view 
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145|| and experimental 
seem to be able 
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Forthcoming satellite missions like PLANCK 
to achieve a level of sensitivity which will enrich decisively our experimental knowledge 
of the CMB polarization with new direct measurements. 

If the background geometry of the universe is homogeneous but not isotropic the 
CMB is naturally polarized [|145|| . This phenomenon occurs, for example, in Bianchi- 



type I models ||148|| . On the other hand if the background geometry is homogeneous and 
isotropic (like in the Friedmann-Robertson- Walker case) it seems very reasonable that 
the CMB acquires a small degree of linear polarization provided the radiation field has a 



non- vanishing quadrupole component at the moment of last scattering [149 . 
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Before decoupling photons, baryons and electrons form a unique fluid which possesses 
only monopole and dipole moments, but not quadrupole. Needless to say, in a homo- 
geneous and isotropic model of FRW type a possible source of linear polarization for 
the CMB becomes efficient only at the decoupling and therefore a small degree of linear 
polarization seems a firmly established theoretical option which will be (hopefully) sub- 
jected to direct tests in the near future. The linear polarization of the CMB is a very 
promising laboratory in order to directly probe the speculated existence of a large scale 
magnetic field (coherent over the horizon size at the decoupling) which might actually 
rotate (through the Faraday effect [||, 0, ^) the polarization plane of the CMB. 

Consider, for instance, a linearly polarized electromagnetic wave of physical frequency 
uj traveling along the x direction in a cold plasma of ions and electrons together with a 
magnetic field {B) oriented along an arbitrary direction ( which might coincide with x in 
the simplest case). If we let the polarization vector at the origin (x = ?/ = 2; = 0, t = 0) 
be directed along the y axis, after the wave has traveled a length Ax, the corresponding 
angular shift (Aa) in the polarization plane will be : 



A« = /e7^f^') (5.x) Ax (8.1) 



2 

'2m \ CO J 

(conventions: ujb = eB/m is the Larmor frequency; Upi = ^J^Tm^e^'/m is the plasma 
frequency rig is the electron density and fe is the ionization fraction ; we use everywhere 
natural units h = c = ks = It is worth mentioning that the previous estimate of 
the Faraday rotation angle Aa holds provided a; 3> a;^ and u ^ Upi. From Eq. (|8.1|) 
by stochastically averaging over all the possible orientations of B and by assuming that 
the last scattering surface is infinitely thin (i.e. that Ax/e^e — Cj.^ where ax is the 
Thompson cross section) we get an expression connecting the RMS of the rotation angle 
to the magnitude of -B at t ~ tdec 



{{Aa)Y'^^-Q'{^^)[—), i?c = 10-^ Gauss, um ^3 xlO'' Hz (8.2) 




(in the previous equation we implicitly assumed that the frequency of the incident electro- 
magnetic radiation is centered around the maximum of the CMB). We can easily argue 
from Eq. ( ^.2|) that if B(tdec)^Bc the expected rotation in the polarization plane of the 
CMB is non negligible. Even if we are not interested, at this level, in a precise estimate 
of Aa, we point out that more refined determinations of the expected Faraday rotation 
signal (for an incident frequency um ~ 30 GHz) were recently carried out ||150| , |151|] 
leading to a result fairly consistent with ( p.l| ). 

Then, the statement is the following, //the CMB is linearly polarized and if a large 
scale magnetic field is present at the decoupling epoch, then the polarization plane of the 



CMB can be rotated ||144|| . The predictions of different models can then be confronted 
with the requirements coming from a possible detection of depolarization of the CMB 
1141. 
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9 Concluding Remarks 



The large scale magnetic fields observed today in the Universe may or may not be primor- 
dial and there could indeed be different possibilities. It could be that in the past history 
of the Universe very strong magnetic fields have been created. These fields could be 
strong enough to affect phase transitions and other phenomena in the life of the Universe 
but, at the same time, too weak to be responsible for the origin of large scale magnetic 
fields. It could also be that magnetic field were indeed strong enough to act as seeds of 
presently observed magnetic fields and, in this case we should be able to find evidence 
that this was indeed the case. In light of this perspective various "observables" , possibly 
affected by the existence of primordial magnetic fields, could be proposed. They include 
the stochastic GW backgrounds, the Faraday rotation of CMB and the baryon asymme- 
try of the Universe. From a more theoretical perspective, primordial magnetic fields can 
be connected to the exsistence of (small and large) extra-dimensions and to the possible 
dynamics of gauge couplings in the early stages of the evolution of the Universe. 
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